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Section  I:  Project  Summary 


1.  Overview  of  Project 

This  project  is  performed  under  the  Office  of  Naval  Research  program  on  Basic  and  Applied  Research  in  Sea- 
Based  Aviation  (ONR  BAA12-SN-0028).  This  project  addresses  the  Sea  Based  Aviation  (SBA)  initiative  in 
Advanced  Handling  Qualities  for  Rotorcraft. 

Landing  a  rotorcraft  on  a  moving  ship  deck  and  under  the  influence  of  the  unsteady  ship  airwake  is  extremely 
challenging.  In  high  sea  states,  gusty  conditions,  and  a  degraded  visual  environment,  workload  during  the 
landing  task  begins  to  approach  the  limits  of  a  human  pilot’s  capability.  It  is  a  similarly  demanding  task  for 
shipboard  launch  and  recovery  of  a  VTOL  UAV.  There  is  a  clear  need  for  additional  levels  of  stability  and 
control  augmentation  and,  ultimately,  fully  autonomous  landing  (possibly  with  manual  pilot  control  as  a  back-up 
mode  for  piloted  flight).  There  is  also  a  clear  need  for  advanced  flight  controls  to  expand  the  operational 
conditions  in  which  safe  landings  for  both  manned  and  unmanned  rotorcraft  can  be  performed.  For  piloted 
rotorcraft,  the  current  piloting  strategies  do  not  even  make  use  of  the  available  couplers  and  autopilot  systems 
during  landing  operations.  One  of  the  reasons  is  that,  as  the  deck  pitches  and  rolls  in  high  sea  states,  the  pilot 
must  maneuver  aggressively  to  perform  a  station-keeping  task  over  the  landing  spot.  The  required  maneuvering 
can  easily  saturate  an  autopilot  that  uses  a  rate  limited  him  system.  For  fly -by-wire  aircraft,  there  is  evidence  that 
the  pilot  would  simply  over -compensate  and  negate  the  effectiveness  of  a  translation  rate  command/position  hold 
control  mode.  In  addition,  the  pilots  can  easily  over -torque  the  rotorcraft,  especially  if  they  attempt  to  match  the 
vertical  motion  of  the  deck. 

This  project  seeks  to  develop  advanced  control  law  frameworks  and  design  methodologies  to  provide 
autonomous  landing  (or,  alternatively,  a  high  level  of  control  augmentation  for  pilot -in-the -loop  landings).  The 
design  framework  will  focus  on  some  of  the  most  critical  components  of  autonomous  landing  control  laws  with 
the  objective  of  improving  safety  and  expanding  the  operational  capability  of  manned  and  unmanned  rotorcraft. 
The  key  components  include  approach  path  planning  that  allows  for  a  maneuvering  ship,  high  performance 
station-keeping  and  gust  rejection  over  a  landing  deck  in  high  winds/sea  states,  and  deck  motion  feedback 
algorithms  to  allow  for  improved  tracking  of  the  desired  landing  position  and  timing  of  final  descent. 

2.  Activities  this  period 

Task  4  -  Dynamic  Inversion  Control  Design 

Work  at  Penn  State  has  been  directed  towards  transitioning  the  Dynamic  Inversion  (DI)  Control  Laws  (CLAWS) 
to  the  Flightlab  Control  System  Graphical  Editor  (CSGE).  PSU  has  developed  scripts  in  Flightlab  to  automate 
generation  of  linear  models  and  synthesis  of  inverse  model  tables.  4th  order  inverse  models  are  generated  from 
trimmed  flight  conditions  starting  at  hover  up  to  160  knots  forward  speed  in  20  kt  increments.  These  are  then 
ported  directly  to  the  CSGE  control  law  diagrams  as  table  look  ups.  Using  these  scripts,  the  re-design  of  the 
control  laws  for  the  various  configurations  (UAV  and  H-53  class)  will  be  relatively  seamless.  In  January  2015, 
the  inner  loop  control  laws  were  completed.  The  inner  loop  CLAWS  include  the  following  response  types: 
Attitude  Command  /  Attitude  Hold  control  laws  in  roll  and  pitch,  yaw  rate  command  /  turn  coordination  control 
laws  in  yaw,  and  vertical  speed  command  control  laws  in  the  vertical  axis.  The  response  types  are  achieved  via 
feed  forward  and  feedback  compensation  schemes  following  the  DI  design  method.  Figure  1  shows  the  top-level 
CSGE  schematic  of  the  inner  loop  DI  controller.  The  control  laws  have  been  tested  and  shown  to  have  good 
tracking  performance  and  good  stability.  Figures  2  and  3,  show  a  response  to  a  yaw  doublet  input  in  hover  and  a 
longitudinal  doublet  input  at  80  knots  forward  speed,  respectively.  The  yaw  doublet  shows  a  180°  heading 


change  in  each  direction.  Off-axis  coupling  is  minimal.  The  aircraft  picks  up  about  2.4  ft/sec  lateral  speed,  but 
this  will  be  corrected  with  the  outer  loop  once  it  is  implemented.  The  pitch  doublet  shows  an  acceleration  to  102 
knots  with  10°  nose  down  attitude,  followed  by  a  deceleration  back  to  80  knots  with  a  10°  nose  up  attitude.  The 
maneuver  demonstrates  the  inversion  model  scheduling  as  the  aircraft  accelerates  across  airspeeds.  There  is 
some  significant  altitude  cross-coupling  with  a  30  ft  altitude  change  with  the  large  pitch  attitude  variations.  This 
could  be  improved  with  a  higher  order  inverse  model  or  through  the  outer  loop  control  law.  Cross-coupling  to 
other  axes  is  minimal  and  the  response  is  very  stable. 


Figure  1  CSGE  Diagram  of  Dynamic  Inversion  Control  Law 


Figure  2  Utility  Helicopter  Simulation  Responses  with  Inner  Loop  DI  CLAW  -  Pedal  Doublet,  Hover 


Nonlinear  time  response 


Inertial  Vy  [ft/sec] 


Figure  3  Utility  Helicopter  Simulation  Responses  with  Inner  Loop  DI  CLAW  -  Long.  Doublet,  80  knots 


Task  5  —  Deck  Motion  Prediction  Algorithm 

During  the  previous  reporting  period,  efforts  were  made  to  implement  a  dynamic  forecasting  method,  a  Holt- 
Winters  (H-W)  method.  In  implementation,  a  multi-block  method  was  used  where  the  data  were  organized  in 
multiple  sampling  blocks  with  offsets.  The  H-W  algorithm  was  then  applied  to  each  data  block  with  adaptation 
to  generate  the  forecasting  for  that  block.  Although  it  provides  a  reasonable  prediction  for  a  short  amount  of 
future  time,  the  algorithm  had  difficulty  predicting  a  longer  time  period  due  to  the  requirement  of  “seasonal 
cycle”  input  which  defines  the  deterministic  trend  term  of  a  given  signal.  Since  the  ship  motion  is  an  inherently 
random  process  (even  though  a  set  of  dominant  periods  can  be  extracted),  the  H-W  algorithm  may  not  be  the 
best  choice  for  the  current  task.  Therefore,  efforts  were  made  in  this  reporting  period  to  investigate  an  alternative 
method,  a  minor  component  analysis  (MCA)  method. 

The  minor  component  is  the  direction  in  which  the  data  has  the  smallest  covariance.  The  statistical  method  for 
extracting  a  minor  component  from  the  input  data  is  called  minor  component  analysis.  The  MCA  determines  the 
directions  of  smallest  variance  in  a  distribution.  They  correspond  to  the  directions  of  those  eigenvectors  of  the 


covariance  matrix  of  the  data  which  have  the  smallest  eigenvalues.  The  main  idea  of  MCA  was  applied  for  a 
curve  fitting  problem  (Oja  1992).  Usually,  the  least  squares  method  is  used  to  solve  such  problems.  For  example, 
given  a  set  of  data  points  (xi,  x2),  the  problem  of  having  a  line  model  to  fit  the  data  in  the  usual  least  square 
sense  becomes  the  problem  of  finding  a  pair  of  estimates.  If  it  is  assumed  that  only  the  measurements  x2  contain 
errors  while  the  measurements  x2  are  accurate,  the  total  least  squares  approach  gives  the  optimal  way  to 
minimize  the  sum  of  the  squared  lengths  of  all  the  bars  which  are  perpendicular  to  the  estimated  line.  The  total 
least  squares  fitting  problem  can  be  reduced  to  the  problem  of  finding  the  minimum  eigenvalue  and  its 
corresponding  normalized  eigenvector  of  matrix  or,  in  other  words,  finding  the  first  minor  component  of  the  data 
set. 

To  implement  the  MCA  method,  a  set  of  ship  motion  data  is  aligned  into  a  sequence  of  vectors,  Xt.  The 

»  1  iV  'T' 

eigenvalues  and  eigenvectors  of  the  autocorrelation  matrix,  R=  [_1  ,  can  then  be  calculated.  The  vector 

Xi  is  formed  as  [Xu,  X2i]T.  where  Xn  is  the  measured  ship  motion  and  X2l  is  the  forecasted  motion  in  the  length  of 
the  forecasting  window.  Based  on  the  MCA  algorithm,  the  forecasted  vector  (X2i)  is  calculated  using  an 
approximated  equality  formulation  consisting  of  eigenvectors  which  are  associated  with  the  smallest  eigenvalues 
of  the  autocorrelation  matrix  (. R ). 

For  simulation  tests  of  the  MCA  based  ship  motion  forecasting  algorithm,  a  set  of  ship  motion  data  is  required. 
Although  it  is  best  to  used  measured  ship  deck  motion  data  for  the  simulation  tests,  the  ship  motion  outputs  from 
USN  SMP  and  STH  are  used  for  the  current  research  because  of  a  lack  of  measured  data.  The  full  6-DOF  motion 
(surge,  sway,  heave,  roll,  pitch,  and  yaw)  is  considered  in  response  to  the  variation  of  the  sea  state  wave 
conditions,  the  wave  heading  angle,  and  the  ship  speed.  Given  a  sea  state,  ship  speed,  and  wave  heading  angle, 
the  ship  motion  is  generated  by  sweeping  the  significant  wave  height  and  the  wave  modal  period  over  the  range 
as  defined  by  a  sea  state  table.  For  this  research,  two  classes  of  ship  model,  DDG-81  class  and  LHA  class,  were 
used  to  generate  the  ship  motion  data.  A  total  of  1,260  test  cases  were  generated  for  each  ship  class.  Each  set  of 
test  conditions  is  a  combination  of  (a)  three  sea  states  (3,  5,  and  6),  (b)  two  ship  speeds  (10  and  20  knots),  (c)  ten 
wave  heading  angles  (0,  ±30,  ±60,  ±90,  ±135,  ±180  degrees),  (d)  five  significant  wave  heights  for  each  given 
sea  state,  and  (e)  four  wave  modal  periods  for  each  given  sea  state.  The  significant  wave  heights  and  the  wave 
modal  periods  were  arbitrarily  selected  from  the  NATO  Sea  State  Numeral  Table  for  the  Open  Ocean  North 
Atlantic  in  order  to  simulate  a  greater  variety  of  ship  motions. 

The  6-DOF  landing  deck  motions  were  sampled  at  a  10  Hz  rate  for  each  test  case.  It  should  be  noted  that 
Landing  Spot  8  was  used  for  the  LHA  class  ship.  The  current  MCA  based  ship  motion  forecasting  algorithm  was 
set  to  predict  the  deck  motion  up  to  6  seconds  ahead  considering  the  accuracy  and  computational  complexity. 
Lurther  research  will  be  required  to  enhance  the  overall  performance.  About  the  first  1,500  points  of  deck  motion 
data  were  used  for  the  MCA  algorithm  training  and  the  remaining  8,500  points  were  used  to  test  the  proposed 
forecasting  algorithm.  The  minor  components  of  the  autocorrelation  matrix  were  selected  in  a  sense  that  their 
energy  added  up  to  no  more  than  1%  of  the  total  energy.  Ligure  4  shows  the  sample  time  history  comparison 
with  6  seconds  of  forecasted  values  of  the  DDG-81  class  ship  for  the  following  cases:  sea  states  3  and  5,  a  ship 
speed  of  20  knots,  and  a  wave  heading  angle  of  0  degrees.  The  black  solid  lines  represent  the  actual  ship  motions 
and  the  red  dashed  lines  are  the  forecasted  values. 

In  order  to  provide  a  quantified  performance  criteria,  several  statistical  terms  were  used.  Lrom  the  initial  testing, 
it  was  observed  that  the  forecasting  error  variations  in  terms  of  standard  deviation  are  very  close  to  a  normal 
distribution.  Thus,  a  range  of  forecasting  error  was  calculated  such  that  the  forecasted  deck  motion  was  within 
90%  of  its  real  value.  Lor  example,  if  a  range  of  forecasted  error  is  wide,  then  the  forecasting  performance  is  not 
acceptable.  In  addition,  the  significant  amplitude,  which  is  based  on  a  concept  similar  to  the  significant  wave 


height,  was  defined  to  represent  the  nominal  level  of  ship  motion.  The  significant  wave  height  is  a  statistical 
term  that  represents  the  average  of  the  highest  1/3  of  the  waves  in  a  given  wave  train.  Since  it  is  known  that 
about  16  percent  of  the  waves  will  be  higher  than  the  significant  wave  height,  it  is  suitable  to  illustrate  the  sea 
condition.  Similarly,  the  ship  motion  can  be  quantified  using  the  significant  amplitude  (e.g.,  significant  heave, 
significant  roll,  etc.).  Figure  5  shows  the  corresponding  quantified  performance  of  the  forecasting  algorithm  for 
the  same  simulation  condition.  Figure  6  shows  the  simulation  results  using  an  LHA  class  ship  for  the  same 
simulation  condition.  It  can  be  seen  that  the  forecasting  error  increases  when  the  sea  wave  heading  angles  are 
close  to  the  port  and  starboard  sides  of  ship  due  to  the  highly  coupled  roll  and  pitch  motion.  This  is  because  of 
the  loosened  correlation  of  the  matrix  R.  In  fact,  one  of  the  default  procedures  of  helicopter  shipboard  landing  is 
to  align  the  ship  with  respect  to  the  wave  heading  angle.  From  the  figures,  the  MCA  forecasting  algorithm 
forecasted  the  ship  motion  reasonably  well  for  wave  heading  angles  between  ±45  degrees,  which  satisfies  the 
requirements  of  the  current  task.  In  addition,  it  can  be  observed  that  the  forecast  algorithm  performs  very  well 
for  relatively  smaller  ship  motions,  such  as  low  sea  state  condition  and  heavier  ship.  It  should  be  noted  that  the 
forecasting  accuracy  is  significantly  improved  for  shorter  forecasting  time  and  it  is  useful  for  a  deck  motion 
feedback  controller  (see  Figure  7). 


Sample  time  responses  -  0  degree  wave  heading 


Sample  time  responses  -  0  degree  wave  heading 


Ship  Type:  DDG-8 1  Class  (Landing  deck)  Sea  state  :  5 

Ship  speed  (knots):  20  Forecasted  time  (sec) :  6 


Ship  Type:  DDG-8 1  Class  (Landing  deck)  Sea  state  :  3 

Ship  speed  (knots):  20  Forecasted  time  (sec) :  6 


(a)  sea  state  :  5 


(b)  sea  state  :3 


Figure  4  Comparisons  of  time  responses  of  forecasted  6-DOF  ship  motion  (DDG-8 1  class) 
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Significant  amplitude  and  MCA  based  forecasted  error  (90%  tolerance) 


Significant  amplitude  and  MCA  based  forecasted  error  (90%  tolerance) 


Ship  Type:  DDG-81  Class  (Landing  deck) 
Ship  speed  (knots):  20 


Sea  state :  5 

Forecasted  time  (sec) :  6 
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Ship  Type:  DDG-81  Class  (Landing  deck)  Sea  state  :  3 


Ship  speed  (knots):  20 


Forecasted  time  (sec) :  6 
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(a)  sea  state  :  5 


(b)  sea  state  :  3 


Figure  5  Quantified  performance  criteria  (vs  wave  heading  angle)  -  DDG-8 1  class 
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Significant  amplitude  and  MCA  based  forecasted  error  (90%  tolerance) 


Signillcant  amplitude  and  MCA  based  forecasted  error  (90%  tolerance) 
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Figure  6  Quantified  performance  criteria  (vs  wave  heading  angle)  -  LHA  class 
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Figure  7  Quantified  performance  criteria  for  3  second  forecast 


Task  6  -  Path  optimization  of  shipboard  helicopter: 

On  the  path  optimization  front,  the  team  is  ready  to  move  forward  with  the  first  path  optimization  study.  This 
study  involves:  (i)  the  mathematical  description  of  a  spatially  and  temporally-varying  approach  profile  using  a 
design  vector,  X,  (ii)  the  development  of  a  novel  objective  function,  F(X),  that  serves  as  a  quantitative 
assessment  of  the  approach  performance,  and  (iii)  mathematical  constraints,  g;(X),  which  are  imposed  to  ensure 
that  the  optimization  results  are  operationally  feasible  and  safe. 

The  details  of  the  formulation  were  described  in  the  previous  progress  report,  and  work  in  the  past  quarter 
focused  on  the  development  of  alternative  objective  functions.  The  overall  performance  for  a  given  approach 
trajectory  can  potentially  be  influenced  by  a  number  of  independent  and  interdependent  performance  factors, 
such  as:  1)  the  maneuver  duration,  2)  the  power  requirements  over  the  course  of  the  maneuver,  3)  path  error 
along  the  maneuver,  and  4)  the  effects  of  the  turbulent  ship  airwake.  Research  to  date  has  focused  on  the 
integration  of  these  factors  into  a  single  objective  function.  While  the  simplest  objective  function  can  be 
constructed  by  applying  equal  weighting  to  all  performance  factors,  alternative  objective  functions  may  be 
constructed  with  varied  relative  weights  between  the  different  performance  factors  (see  Figures  below).  The 
representative  results  [from  our  abstract]  shown  here  are  based  upon  a  simple  sensitivity  study  that  will  be 
extended  for  a  full  optimization  in  the  coming  quarter. 


Figure  8  Sensitivity  Study  of  a  Weighted  Performance  Index  for  Automated  Ship  Approaches  (These  contour 
plots  show  the  proposed  weighted  cost  function  for  the  automated  ship  approaches  presented  in  the  previous 
progress  report,  the  cost  is  a  function  of  the  approach  azimuth,  glide  slope,  initial  velocity,  and  aggressiveness  of 

deceleration). 


3.  Significance  of  Results 

The  deck  motion  forecast  method  was  established  using  the  MCA  algorithm  and  was  tested  for  various 
conditions  which  were  combinations  of  multiple  sea  states,  ship  speed,  wave  heading  angle,  significant  wave 
height,  and  average  wave  modal  period.  The  simulation  results  identified  the  capability  of  the  MCA  based  ship 
motion  forecasting  algorithm  by  providing  a  forecast  confidence  measurement  in  terms  of  a  statistical 
interpretation  of  the  prediction  error. 

Implementation  of  the  DI  Controller  in  CSGE  /  Flightlab  has  so  far  been  successful.  The  design  scripts  are 
automated  so  that  implementation  for  UAV  and  H-53  class  rotorcraft  can  be  performed  very  quickly  (just  a  few 
hours  work  to  integrate  and  run  inverse  model  generation  for  preliminary  implementation,  followed  by  a  week  of 
work  for  testing  and  tuning). 

4.  Plans  and  upcoming  events  for  next  reporting  period 


Task  1  -  Plant  Model  and  Disturbance  Models'.  The  FireScout  and  H-53  class  FLIGHTLAB  models  will  be 
developed  and  distributed  to  team  members. 


Task  4  -  Dynamic  Inversion  Control  Design'.  Will  complete  implementation  and  testing  of  outer  loop  DI 


CLAWS.  A  preliminary  version  of  the  CLAWS  will  be  transitioned  to  counterparts  at  NAVAIR  and  NSWCCD, 

Task  5  -  Ship  Motion  Prediction'.  Next  effort  will  be  focused  on  the  MCA  algorithm  enhancement  to  improve  the 
accuracy. 

Task  6  -  Path  optimization  of  shipboard  helicopter:  Continue  with  a  full  path  optimization  study  following  the 
plan  in  our  ALAA  draft  manuscript.  We  are  also  likely  to  continue  path  optimization  work  using  our  Flightlab 
simulation  once  it  is  complete  (as  opposed  to  the  GENHEL  simulations  used  to  date). 
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